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We study the survival probability time distribution (SPTD) in dielectric cavities. In a circular 
dielectric cavity the SPTD has an algebraic long time behavior, ^ t~^ in both TM and TE cases, but 
shows different short time behaviors due to the existence of the Brewster angle in TE case where the 
short time behavior is exponential. The SPTD for a stadium-shaped cavity decays exponentially, 
and the exponent shows a relation of 7 ^ 1 is the refractive index, and the proportional 
coefficient is obtained from a simple model of the steady probability distribution. We also discuss 
about the SPTD for a quadrupolar deformed cavity and show that the long time behavior can be 
algebraic or exponential depending on the location of islands. 
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I. INTRODUCTION 

Recently, lasing modes from dielectric microcavities 
have attracted much attention due to its potential ap- 
plication to optoelectric circuits and optical communica- 
tions In particular, there was a lot of theoretical and 
experimental effort to excite directional lasing modes in 
deformed microcavities |^J. It is now well known that 
the lasing pattern has a very close relationship with the 
internal ray dynamics given by the boundary geometry 
of cavity. It is also reported that the property of the 
openness of the dielectric cavity plays an important role 
in the resonance pattern analysis (J, U] . 

For a general open system, the survival probability 
time distribution (SPTD) or its time derivative, the es- 
cape time distribution, is a basic physical quantity. Many 
studies are focused on the relation between the long time 
behavior of the SPTD and the internal dynamics, and 
it is known that the SPTD has algebraic and exponen- 
tial decays in integrable and chaotic systems, respectively 
H, 0, 0, M S mixed systems, having both inte- 

grable islands and chaotic sea in phase space, the SPTD 
has algebraic long time behavior which originates from 
the slow escape mechanism due to the stickiness of K AM 

tori ini. 

The property of openness of the dielectric cavity 
is different from the open systems previously studied 
5, 6, 7, ^05^3- Rays can escape through any boundary 
point, and partial escapes, depending on the incident an- 
gle, are possible. This unique property can be reflected 
on the long time behavior of the SPTD. 

In this paper, we study the SPTD in dielectric cavities 
of various boundary geometries such as circle, stadium, 
and quadrupole, which are typical examples of integrable, 
chaotic, and mixed systems, respectively. The SPTDs in 
these dielectric cavities show basically similar behavior to 
the open cavity with a small hole on the cavity bound- 
ary, but the exponents are different. In particular, we 
show that the ergodic property cannot be applied for the 
stadium-shaped dielectric cavity even in the small open- 
ing limit, n ^ 00, n is the refractive index. 



The paper is organized as follows. In Sec. II the al- 
gebraic long time behavior of the SPTD in the circu- 
lar dielectric cavity is derived analytically and confirmed 
numerically for both TM and TE waves. It is shown 
in Sec. Ill that the SPTD for a stadium-shaped cavity 
decays exponentially, and the exponent 7 has ~ de- 
pendence and the proportional coefficient can be under- 
stood from a simple model of the steady probability dis- 
tribution (SPD). The SPTD in the quadrupole-deformed 
dielectric cavity is discussed in Sec. IV and we finally 
summarize the results in Sec. V. 



II. CIRCULAR DIELECTRIC CAVITY - 
INTEGRABLE SYSTEM 

Many authors have studied the SPTD for the open 
billiard with a small hole on boundary H, IE S 0- It 
is known that the SPTD in a circular billiard decays al- 
gebraically, Psv{t) oc t~^. In this section we study the 
SPTD for the circular dielectric cavity, and it will be 
shown that the SPTD shows still an algebraic decay but 
the exponent is different. 

For simplicity, we focus on TM wave first and for TE 
wave we will mention only difference later. In the circular 
geometry, ray dynamics is integrable and rays in the open 
area of the phase space, i.e., —l/n<p< 1/n, p = sin^, 
6 being the incident angle, can partially escape from the 
cavity. The formal expression of the SPTD in the circular 
dielectric cavity is given by 



Psv{t) = 



2L 



ds / dp R{p) 



Nit) 



(1) 



where L is the boundary length, and Pc is the critical line 
for total internal reflection, i.e., Pc = sin^c = l/n, and 
R{p) is the reflection coefficient for TM wave [13, 



R{0) 



n cos — cos Ot 
n cos + cos 6t 



(2) 
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where nsin^ = sin^^, and N(t) is the number of bounce 
on the boundary. Since N{t) = t/2cos0 in the circular 
geometry, when considering a unit circle and a time scale 
as the length of ray trajectory, Eq.Q can be rewritten 



where 



Psv{t) 



G{0) 



■ n 
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Jo 



dO cos(9exp[-G((9)t], 



cos 



In 1 



2 cos Of 



n cos — cos Ot 



(3) 



(4) 



Note that the rays near the critical line Pc can survive 
longer time and dominate long time tail behavior. There- 
fore, we can expand G{0) from Oc by changing variable, 
6> = 6>c - X, as 



G{e) ^ ax 



1/2 . 



3/2 



where 



and 



a 



2nV2V^ 



(3 = -a{ 



n 



1 



(5) 



(6) 



(7) 



Substituting the lowest term in Eq.(|5|) into Eq.(|3|), we 

can obtain the long time behavior of the SPTD as 



Psv{t) 



2Vn2 



ct)e~ 



2^fr. 



(8) 



We emphasize that the SPTD for the circular dielectric 
cavity decays as as shown in Eq.Q, different from 
the open billiard with a small hole where decays as . 
This means that the property of openness can change the 
exponent of the algebraic decaying SPTD. 

For TE wave case the reflection coefficient is given by 



Rte{0) 



n cos 6t — cos 6 
n cos 9t + cos 



(9) 



and the expansion of G{0) and the SPTD at a long time 
are the same as Eq.© and Eq.Q with different expan- 
sion coefficients, i.e.. 



O-TE 



2nV2Vr?^ 



and 



CtTE 



'TE 



A\fn? 



1 



6). 



(10) 



(11) 



We note that the dependence of the SPTD on the refrac- 
tive index n in TM and TE waves is quite different, i.e.. 




FIG. 1: (color online) The SPTDs of the circular dielectric 
cavity for (a) TM and (b) TE waves. Black circle, red square, 
green diamond, blue triangle (up), and brown triangle (down) 
are for n=2, 4, 6, 8, and 10, respectively. These shows 
long time behavior and a very good agreement with the solid 
lines which represent the n — 2 case in Eq.(|H|. 



Psv{t) ^ nH-'^ for TM case, but Psv{t) ^ n-'^t''^ for 
TE case. The proportionality of of the TE case does 
not mean that the circular cavity with a higher n is more 
leaky, since we take into account of only the open region 
in the phase space, — 1/n < p < 1/n (see Eq.Q). 

In order to perform numerical calculation for the 
SPTD in the circular cavity, we take 10^ random ini- 
tial points in the open region of the phase space. We 
then trace each point with a weight determined by R{p) 
when bouncing from the boundary, and sum the weights 
between t and t + A, we take A = 1 in the calculations, 
for all points in the ensemble, and finally normalize to be 
unit when t = 0. Figure 1 shows the numerical results of 
the SPTD in the circular cavity for TM and TE cases. It 
is clear that the SPTD for both cases shows an algebraic 
long time behavior, ~ and the dependence on n is 
correctly described by Eq.Q which is indicated by the 
solid lines for n = 2 in Fig.l. 

A substantial difference between the TM and TE cases 
appears in the short time behavior. As shown in Fig. 1 
(a) the short time behavior of the TM SPTD is smoothly 
connected to the t~'^ long time tail, on the other hand 
that of the TE SPTD shows rather an abrupt transi- 
tion to the algebraic long time tail and the detail of the 
short time behavior seems to be characterized by an ex- 
ponential decay. This exponential short time decay is 
clear in Fig. 2 (a) and the exponent 7(n), when fitted 
as exp(— 7(n)t), appears as the solid dots in Fig. 2 (b). 
These exponents are well described by a simple approxi- 
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FIG. 2: (color online) (a) The exponential short time behavior 
of the SPTDs for TE case. The colors for different n is the 
same as in Fig. 1. These can be fitted as exp(— 7(n)t) in 
a short time range, (b) The exponents j{n). The solid line 
represents the result of a simple approximation of Ea. (|13|) . 



mation for the reflection coefficient Rte{0). If we expand 
Rte{^) at ^ = and take only the lowest term, then 

Psy{t) - exp(- In(^^)t) = exp(-7(n)t). (12) 
n — 1 

The solid line in Fig. 2 (b) represents the relation 

1. 



7(n) = ln( 



n - 



1' 



(13) 



Here we emphasize that although the lowest term of the 
expansion of Ea. (|T2j) is the same for both TM and TE 
cases, only TE case allows the exponential short time 
decay. The reason for this is the existence of the Brew- 
ster angle in the TE case, Ob = arctan(l/n) where rays 
can escape without reflection, i.e., R{Ob) = 0. The rays 
with the incident angle in the range of —Ob < < Ob 
dominate the short time behavior, while the other parts, 
Ob < \0\ < Oc^ mainly contribute to the long time alge- 
braic tail. 

It is important to know when the algebraic decay starts 
to appear in both TM and TE cases. In TM case, the 
main factor for the deviation from the t~'^ decay comes 
from the finite integral bound, and it corresponds to the 
terms containing the upper bound Oc in Eq.((H|). We then 
estimate the transition time when the deviation from the 
decay is about 10 %, and the result is 



1.38(n^ - 1)^/^ 
nY^arcsin(l/n) 



(X n. 



(14) 



In Fig. 1 (a) the corresponding transition times are in- 
dicated by arrows and show a good agreement with the 



numerical calculations for various n. 

Due to the existence of the Brewster angle, the transi- 
tion time for TE SPTD can be determined by a different 
way. As mentioned above, the TE SPTD shows a short 
time exponential behavior and a long time algebraic be- 
havior. Therefore we can estimate the transition time by 
finding the intersection time for both different behaviors. 
From Eq.(|H|) with arE and Ea. (|T2j) . for a large n we can 
get an implicit equation for the transition time as 



n n 2n^ 



(15) 



The transition times for various refractive indices are in- 
dicated by arrows in Fig. 1 (b) and well represent the 
transition times of the numerical results. The solution tc 
of the above equation cannot be described by a simple 
power of n, but we can show 

tc(n) (xn^(^\ /i(n) > 1. (16) 

If we take a logarithm of Eq. (|T^ , then we get 



tc 

n 



In — = In 2n , 

n 



(17) 



which generally has two solutions and the larger solution 
is relevant. The point to, at which the slopes of the two 
functions in the left hand side of the above equation are 
identical, should locate between the two solutions. By 
differentiating the above equation, we get to = n. There- 
fore, 



tc> to = n. 



(18) 



Even though both TM and TE cases show the same 
t~'^ long time decay in the circular cavity, the short time 
behavior and the n dependence of the transition time 
are quite different. We emphasize that these differences 
originate from the existence of the Brewster angle in the 
TE case. 



III. STADIUM-SHAPED DIELECTRIC CAVITY 
- CHAOTIC SYSTEM 

As an example of chaotic dielectric cavities, we take 
a stadium-shaped one with parallel linear segments of a 
length / and two semicircles of a radius R. The stadium- 
shaped billiard has been a typical chaotic system in the 
research of classical and quantum chaos. The escape 
property through a small hole on the boundary of the 
stadium-shaped billiard has been investigated by many 
authors 9] . They have shown that the escape time distri- 
bution exponentially decays first and later becomes alge- 
braic, and the transition time tc increases as the hole size 
decreases. The algebraic decay at long times comes from 
the stickiness near the marginally stable line in phase 
space corresponding to the bouncing ball trajectories. On 
the other hand, in the stadium-shaped dielectric cavity 
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FIG. 3: The steady probability distributions, (a) n = 2 case 
with system parameters {R,l) — (1,2) for TM wave, (b) 
n = 10 case with system parameters (i?, /) = (1,2) for TM 
wave. The black points indicate the rays traced from a ran- 
dom ensemble with weights higher than 0.1 at about t — 12 in 
n = 2 case and t = 37 in n = 10 case. The solid lines denote 
the critical line =bpc for total internal reflection. 



the ray trajectories of the bouncing ball type cannot con- 
tribute to the long time behavior due to the property of 
openness, i.e., rays with almost vertical incidence escape 
easily and contribute to the short time behavior. As a 
result the SPTD shows only exponential decay (see Fig. 
4, 5). 

The exponential decay in the dielectric chaotic cavity 
implies the existence of the steady probability distribution 
(SPD), Ps{s^p) which is defined as the spatial part of the 
survival probability distribution (5, p, t) 3]. With this 
SPD, we can express the SPTD as 



Psv{t) 



ds i dp Psv{s,p,t) ^ Cexp(-7t), (19) 
J -I 



where C is a constant and 



\ ds I dp Ps{s,p)T{p), 
J-i 



(20) 



where the transmission coefficient T{p) is given as T{p) 
l-R{p). 




FIG. 4: (color online) The SPTD and the decay rates for 
TM wave, (a) The exponential SPTD for a stadium dielectric 
cavity with (i?, /) = (1, 1). Black, red, green, blue, and brown 
lines are for n=2, 4, 6, 8, and 10, respectively, (b) The decay 
rates with increasing n. Black circle and red square represent 
results for system parameters (i?, /) — (1,1) and (1,2), re- 
spectively. The solid line shows behavior. The coefficient 
D is given as D ■ 



Note that the above equation is satisfied in the expo- 
nential decay region and cannot describe the nonexpo- 
nential very short time behavior. From Eq.(|2n|), if we 
know the SPD, we can estimate the decay rate 7. How- 
ever, the structure of the SPD is usually very compli- 
cated because it depends on the openness as well as the 
boundary geometry of the cavity. Figure 3 (a) shows the 
approximate of the SPD when n = 2 which is a snap shot 
of the Psv{s^p^t) captured at about t = 12. The partial 
escape property of the dielectric cavity allows for rays 
to distribute on unstable manifold structure in the open 
region, —l/n<p<l/n. 

Even though it is difficult to estimate the SPD in usual 
cases, for the large n case, the small opening case, we can 
simplify the SPD by assuming a uniform distribution over 
the whole phase space except the open regions related 
to the linear segments of the stadium boundary. The 
approximate of the SPD for n = 10 shown in Fig. 3 
(b) supports this assumption. We note that this is a 
substantial difference from the escape through a small 
hole on boundary where entirely uniform distribution is 
assumed due to the ergodic property 5]. Based on the 
assumption of the partial ergodicity, we can rewritten the 
decay rate as 



7 



ttR 



y/A{L-2l/n 



I dpT{p), (21) 



where we insert the factor 1/a/A, A = ttR^ + 2RI be- 
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FIG. 5: (color online) The SPTD and the decay rates for TE 
wave. Detail of caption is the same as Fig. 4 except D — 




t 

FIG. 6: (color online) The SPTDs in the quadrupolar de- 
formed dielectric cavity. The SPTDs are algebraic in the 
s = 0.1 cases while exponential in the e — 0.2 cases. The 
black and red lines denote the case of TM and TE waves, 
respectively. 



IV. QUADRUPOLE-DEFORMED DIELECTRIC 
CAVITY - MIXED SYSTEM 



ing the area of the stadium, from the consideration of 
time scale. The integral of the above equation means the 
degree of openness and in the large n limit decreases as 
- 2z/7rn-2 for both TM {u = 1) and TE {u = 2) cases 
(see Appendix). Therefore, for the large n limit the decay 
rate becomes 



7 ^ — 1= — n 

Val 



(22) 



Numerical results for the SPTD in the chaotic stadium- 
shaped dielectric cavity are shown in Figs. 4 and 5 for 
TM and TE cases, respectively. We take two systems; 
one is (R^l) = (1,1) and the other is (R^l) = (1,2). 
In calculation, we use a random ensemble of 10^ initial 
points spread over the whole phase space and trace the 
survival probability with time, the time is scaled as the 
length of trajectory in the spatial space as before. The 
exponential behavior of the SPTD is clear even at long 
time limit in both TM and TE cases. This means that 
the sticky region locating on the center of the open region 
in phase space dose not contribute long time decay due 
to its easy escape. The dependence of the decay rate 7 
on the refractive index n shows a very good agreement 
with Eq.(|22) for large refractive indices, in both systems 
with different area A. This implies that even in the small 
opening limit, n ^ 00, we cannot use the ergodic prop- 
erty over the whole phase space. Instead, we have to 
consider structure of the SPD even in the small opening 
limit. 



The escape property in generic mixed systems, show- 
ing a mixed phase space portrait: integrable islands in 
a chaotic sea, has been extensively studied. It is well 
known that the long time behavior of the SPTD is alge- 
braic due to the stickiness of the KAM tori surrounding 
islands 111. 



(23) 



However, there is no rigorous theory expecting the value 
of the exponent 77 which has been estimated based on 
numerical calculations and seems to be nonuniversal. 

In this section, we consider a quadrupolar dielectric 
cavity which is the typical example of a deformed mi- 
crocavity and shows a mixed dynamics. The boundary 
equation is, in the polar coordinates. 



r{(j)) = 1 + £ cos 2^, 



(24) 



where e is the deformation parameter. Here, we present 
numerical results of the SPTD and show that the long 
time behavior of the SPTD is determined by whether 
islands locate in the closed region, pc < \p\ < or not. 

For n = 2 case, we numerically calculate the SPTDs at 
two deformation parameter values, e = 0.1 and 0.2, which 
are shown in Fig. 6. In the e = 0.1 case, the SPTD shows 
an algebraic decay, i.e., Psv{t) ~ which is consistent 

with the previous studies on mixed systems. However, in 
the s = 0.2 case, the long time behavior of the SPTD is 
exponential, i.e., Psv{t) ~ exp(— 0.05t). This clear dif- 
ference of the SPTD between s = 0.1 and 0.2 cases can 
be explained by the phase space portraits. Figure 7 (a) 
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very diflFerent short time behavior due to the existence of 
the Brewster angle in TE case where exponential short 
time behavior is shown. The SPTD for a stadium-shaped 
cavity decays exponentially, and the exponent has a close 
relation to the steady probability distribution (SPD). In 
the large n limit, the SPD can be approximated by an 
assumption of a partial ergodicity, a uniform distribu- 
tion over a specific part of phase space, which gives a 
correct description of the exponent in both TM and TE 
cases. We have also discussed about the SPTD for the 
quadrupolar deformed cavity and shown that the long 
time behavior can be algebraic or exponential, depend- 
ing on the location of islands. 
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FIG. 7: The phase space portraits for the quadrupolar de- 
formed cavity, (a) The s = 0.1 case, (b) The s = 0.2 case. 
The dashed lines denote the critical line pc — 1/2 for total 
internal reflection for n — 2. 



shows the phase space portrait for e = 0.1 case. There 
are many islands in the closed region, < P < 1, so 
the stickiness of the KAM tori delays the ray escape and 
results in the algebraic tail. On the other hand, as shown 
in Fig. 7 (b) there is no island in the closed region for 
£ = 0.2 case and all islands exist in the open region. The 
rays trapped by the stickiness of the KAM tori contribute 
to the short time escape behavior and the resulting SPTD 
shows exponential long time decay. Therefore, the posi- 
tion of islands plays important role to understand the 
SPTD of mixed systems. 

V. SUMMARY 

We have investigated the survival probability time 
distribution (SPTD) in three dielectric cavities show- 
ing different ray dynamics; circle (integrable), stadium 
(chaotic), quadrupole (mixed) shapes. In the circular 
dielectric cavity the SPTD has an algebraic long time 
behavior, ~ in both TM and TE cases, but shows 



Appendix 

Here we present the analytical expression of the degree 
of openness (see Eq.(|23)) for TM wave. The degree of 
openness is defined by 

['^^ dpT{p), (25) 

J-l/n 

where T{p) = 1 — R{p)^ R{p) is given in Eq.(|2|). This 
integral can be expressed by an analytical function as 

(n2 - 1)2^ 4' 2^ ^ 2' 2' ' 

r.A 5,^, 1 1 1 , 40 8 1, 
+i^(2>2)^(-2'2'3;^)--n+--] 

27^n-^ (26) 

where B{x^ y) is the beta function and F(a, 7; z) the 
Gauss hypergeometric function 13]. 

For TE wave, only difference is the replacement of R(jp) 
by Rte{p) of Eq.Q, and the result is 

/ 47rn-2 (27) 

for the large n limit based on a numerical calculation. 
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